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We display a logarithmic divergence in the density matrix of a scalar field in the presence 
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particle. 

" Theory Division, CERN, CH-1211, Geneva, Switzerland, 

^ University of Oxford, Dept. of Physics (Theoretical Physics), 1 Keble Road, Oxford 
0X1 3NP, United Kingdom, 

Center for Theoretical Physics, Dept. of Physics, Texas A & M University, College 
Station, TX 77843-4242, USA and Astroparticle Physics Group, Houston Advanced 
Research Center (HARC), The Mitchell Campus, Woodlands, TX 77381, USA. 
* P.P.A.R.C. Advanced Fellow. 



CERN-TH/95-336 
CTP-TAMU-04/95 
ACT-02/95 
OUTP-96-05P 



January 1996 



1 Introduction and Summary 



There is currently much debate whether microscopic black holes induce quantum 
decoherence at a microscopic level. In particular, it has been suggested [Q that 
Planck-scale black holes and other topological fluctuations in the space-time back- 
ground cause a breakdown of the conventional S'-matrix description of asymptotic 
particle scattering in local quantum field theory, which should be replaced by a 
non-factorizable superscattering operator ^ relating initial- and final-state density 
matrices: 



It has further been pointed out that, if this suggestion is correct, there must be a 
modification of the usual quantum-mechanical time evolution of the wave function, 
taking the form of a modified Liouville equation for the density matrix [Q]: 



The extra term in (0) is of the form generally encountered in the description of 
an open quantum-mechanical system, in which observable degrees of freedom are 
coupled to unobservable components which are effectively integrated over, which 
may evolve from a pure state to a mixed state with a corresponding increase in 
entropy. The necessity of a mixed-state description is generally accepted in the 
presence of a macroscopic black hole, but is far from being universally accepted in 
the case of microscopic virtual black-hole fluctuations. 

Three of us have been analyzing the possibility of quantum decoherence in a non- 
critical formulation of string theory |Q, in which the two-dimensional string black 
hole solution is used as an example of a quantum fluctuation in space time. We 
exhibited in this model extra logarithmic singularities, associated with transitions 
between different conformal field theories on the world sheet, which induced an ex- 
tra term in the quantum Liouville equation of the form conjectured in (^. This 
derivation of @) was, however, incomplete, in that two dimensions are not the same 
as four, and even in two dimensions one does not have a complete non-perturbative 
formulation of string theory that includes transitions between different conformal 
field theories^ In parallel, two of us have been studying a scalar field in the pres- 
ence of a four- dimensional Einstein- Yang-Mills black- hole background [0], and have 
demonstrated that its entropy exhibits a logarithmic divergence that cannot be ab- 
sorbed into a simple renormalization of the Bekenstein entropy of the black hole. 
This can be understood as a quantum reflection of the entanglement of the external 
scalar field with internal black hole properties that alter with its radius. 

The purpose of this paper is to show that this logarithmic divergence, which also 
appears in the partition function Z but not at the same order in the Hamiltonian, 

"'^Progress in this direction is being made via studies of duality and conifolds with consider- 
able similarity to the approach taken in p : see . 





dtp = i[p,H]+^Hp 



(2) 



induces an extra non-commutator term into the quantum Liouville equation for 
the density matrix of the scalar field. This follows from the identification of the 
renormailization scale with time, as in the two-dimensional non-critical string 
and constitutes the first example of a derivation of equation (Q) in four dimensions. 
We argue that, once one allows for the quantum decay of the black hole, the loga- 
rithmic divergence found in the previous static case becomes a new source of time 
dependence for the scalar field that leads to a monotonic increase in entropy, and 
hence cannot be absorbed within the conventional Liouville equation. The order of 
magnitude that we find for the extra term is 

where n is the mass of a quantum of the scalar field. The estimate (0) is the 
maximum that had been indicated by previous string black hole studies 0, and is 
not very far from the sensitivity of present experimental searches M. 



2 Entropy, Hair and Logarithmic Divergences 

As a first step in this analysis, we now review relevant previous work and recall 
some of the ingredients we use. The observation that quantum black holes must 
be described by mixed quantum states goes back to the work of Bekenstein and 
Hawking [|1^, who showed at the tree level that a generic black hole has non-zero 
entropy related to the area of its event horizon: 

This entropy represents the number of quantum black hole states that are not dis- 
tinguishable by traditional measurements using the varieties of hair known in con- 
ventional local four-dimensional quantum field theories. It is known that string 
theories possess an infinite number of local symmetries, and we have suggested that 
a corresponding infinite set of conserved quantum numbers, baptized ly-hair ||TT[| , 
characterize string black hole states, at least in two dimensions. We have also pointed 



out [|T2[ that the number of string black hole states grows in the same way as the 
logarithm of the Hawking-Bekenstein entropy and suggested ways in which the 
W quantum numbers could in principle be measured and used to distinguish these 



string black hole states |11 



However, we have also pointed out that a complete set of these measurements is 
not possible in practice, and have gone on to suggest that the incompletenes of the 
Pi^-hair measurements provides a measure of the loss of information associated with 
a string black hole [0]. The observable low-energy degrees of freedom are entangled 
by W symmetry with unobservable internal states of the black hole. Integration over 
these unobserved degrees of freedom provides a truncated subtheory which possesses 



logarithmic divergences absent in a conformal field theory. The latter would repre- 
sent string theory in a fixed classical background, in which scattering is described 
by a conventional S matrix. The extra logarithmic divergences are associated with 
fiuctuations in the black-hole background, which require transitions between dif- 
ferent conformal field theories. Two-dimensional examples of of these divergences 
are provided by string world-sheet monopoles and instantons, which represent the 
creation and mass shifts of black holes, respectively |Q . After identification of the 
renormalization scale appearing in the logarithms with a time-like Liouville field, 
and thence with the target time variable, we have been able to derive (0) with an 
explicit representation for the non-commutator term. 



We now point out the analogies with the study of a scalar field coupled to an 
Einstein- Yang-Mills (EYM) black hole in four dimensions. This study was purely 
in the context of quantum gravity, with no string degrees of freedom. On the other 
hand, the presence of the gauge field endows the EYM black hole with additional 
gauge hair (c.f., the W hair of the two-dimensional string black hole) that is entan- 
gled with the scalar field via gauge interactions (c.f., non-diagonal W generators). 
The renormalizability of the EYM field theory enables the partition function and 
entropy of the four- dimensional black hole to be calculated with the inclusion of 
quantum corrections, with the results 
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where Sbh is the tree- level Bekenstein-Hawking entropy (^), and the Sq are the 
quantum corrections, indicated by square brackets. The subscript h denotes quanti- 
ties at the horizon of the black hole, is the horizon radius, the symbols X denote 
ratios X/r/j, primes denote differentiation with respect to the radial coordinate r, 
m is the mass function, 5 is defined in ref. and will be given explicitly below, 
and e is a small, positive fixed distance which will play the role of an ultra-violet 



cut-off. Its presence is associated with the 'brick wall' boundary condition |T3| for 
the wavefunction of the (scalar) matter fields in the black hole background. 



The first term in the quantum corrections to the entropy (|^) can be absorbed 
into the bare lowest-order Hawking-Bekenstein entropy (|^) via a renormalization 
of Newton's constant Gn : G^^ Gjy^ren- However, the second quantum correc- 
tion cannot simply be absorbed into a bare parameter in this way, and therefore 
corresponds to a new effect beyond the reach of the conventional renormalization 
programme, as we discuss in the next section. Formally, it may be absorbed into 
a "quantum" horizon area, corresponding to a "quantum" radius r^, as we shall 
discuss in section 4. 



3 Hamiltonian of the Scalar Field in an EYM 
Black-Hole Background 

In order to gather information needed for the interpretation of the logarithmic 
divergence in the partition function @ and in the entropy @ 0, we now compute 
the Hamiltonian of a scalar field of mass fi described by the Lagrangian 

= -lV^{9'%<pdx<P + fi'<P') (7) 

which satisfies the Klein-Gordon equation in an EYM black-hole background de- 
scribed by the metric 

ds^ = -e^dt^ + e^dr^ + r\de^ + sin Od^^) (8) 

where the forms of the metric functions r,A are discussed in []^ and below. It is 
easy to verify that the canonical momentum vr which is conjugate to (p is given in 
this case by 

^ = 7^ = -v^^?*'^^^'^ = e^(-^+^V2sin^9,0 (9) 
d[dt(p) 

Defining the Hamiltonian density = vr^u^ — £ in the normal way, we arrive at 
the following expression for the Hamiltonian H: 

H = J drdOdipH : 

\{de(j)f + \cosec'e{d^(j)f + fi^(j)') (10) 



which we now evaluate. 



To evaluate (0), we first expand <^ in normal modes 

0(t, r, ^) = e-''''fEi{r)ZEi{9, ^) (11) 



where Zim is a real spherical harmonic 



, and fEi satisfies the radial equation 



dr 



e5(r-A)^2 



d 



r —Jei 
dr 



2 /i?; = (12) 



which we solve with the 'brick-wall' boundary conditions |]T3| = at r = r/i + e 
and r = L : L » Th- Here r/j is the horizon radius of the black hole, e is a small, 
positive fixed distance which will play the role of an ultra-violet cut-off, and L is 
an infra-red cut-off. Introducing annihilation and creation operators a and for 
each normal mode (pJ]) in the usual way, we arrive at the following expression for 
the Hamiltonian: 



1 



-E + E^Jei 



(13) 



E,l,m 

where all the dependence on e is contained in the last factor in (^) and 
Before we estimate this integral, a few comments are in order concerning ([T3|). 



Jei = r rV^^^'-'^^^fl^dr. 

Jru+e 



(14) 



First, the form of the Hamiltonian ensures that the Fock basis states are eigen- 
states of energy, so that the Hamiltonian can be written in the alternative form 



states 

X I ^Uk, , ^nfc2 , . . . . ) {^Uk^ , ^nfca , • • • 



(15) 



where the expectation values are, from (p!3|). 



modes 
2 4" 



1 



-E + E'J, 



El 



^E, + E 



Eih 



(16) 



The first term in (|T6[) is state-independent and contributes an infinite constant to 
the partition function Z = Tie"^^ . We shall ignore this factor from here on, since 
it does not contribute to correlation functions. The second term in (^) will be 
compared later in this section to that used in to calculate the partition function. 

In order to estimate the integral Jei, we shall use the WKB approximation for 
the /ei: 



fEi{r) 



1 -ifr-Ai ^Ei 
-e 2^ ' , -sm 



v^(r) 



KEi{r')dr' 



(17) 



where the radial wave number K is given by 



;i8) 



The constant Aei in each case is determined from the normahzation conditions for 
the /'s, namely 



and the boundary condition 
implies that 



fEi{L) = 



/ KEi{r')dr' = rni for some integer n. 



Next introduce a dimensionless co-ordinate x given by 

r 



X 



and the rescaled function K as follows: 

Kei = thKei 

and 

e 



(19) 
(20) 
(21) 

(22) 

(23) 
(24) 



We decompose the integral in ([19D as A\i{Ii + 12), where /i is the contribution to 
the integral for x very close to one, and I2 is the remainder of the integral, so that 

^ (25) 



42 



2E{h + h) 

Similarly, we may write Jei = A'^iiJiEi + J2Ei)- 



Writing the metric functions in the form: 



2m{r) 



-25(r) 



.A 
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2m{r)\ ^ 



r 



it can easily be seen that, for spacetimes in which 

e-^^^*-) = 1, 



(26) 



(27) 



the integrals (|Tj) and (0) are identical. This is the case for a black hole without 
hair, i.e., the Schwarzschild solution. In this case, Jei = '^/2E, the Hamiltonian H 
(p!3D is independent of loge, and the form of the partition function Z is identical to 
the expression (|^) found in M. 



For geometries with hair, the contribution to the integral in ( [T^ ) for x ~ 1 is 
dependent on log e, and is given by 



J 



lEl 



^El 
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(28) 



For other values of x, the integrand contains a factor of the generic form 



2ruE 
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(30) 



which is the same for both I2 and J2. Note that the logarithmic terms in the 
expressions ( p^ ) and ( pO] ) have as their argument e = e/vh rather than L/e as has 
been discussed by |jl5|, the contribution to the entropy for large values of r being 
proportional to The reason for this is that the expansion we are using for 



the integrand in (|1^) is valid only for values of x very close to 1. An alternative 
expansion in inverse powers of x has to be used for large values of x. 



We therefore write 

I2 = CT)sm^ 
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(31) 



where the *jF's and ^Q^s are constants which are independent of e but do depend on 
L. In total, therefore, we have 



J, 



1 J1EI + J2EI 



El 



2Ee^^^JiEi + hEi 
1 



2E 



2S^ . J-2EI- hEie 
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e^^^JiEi + hEi 



(32) 



The second term in this expression is of the order of (loge) ^, and hence is very 
small. 



We now proceed to evaluate the rate of change of Jei with respect to loge. The 
result is: 

dJEi 11" 



a loge 2Ee''^^JiEi + I2EI 



dJ2El _ ^-2&h ^^"^El 



+ smaller terms. (33) 



d log e d log e 

We therefore arrive at the following expression for the e-dependence of the scalar- 
field Hamiltonian H: 

dH 1 ^-^ . + + . „o dJf 



To estimate the order of magnitude of (0), we need to estimate the coefficients 
J-' and Q. It is not possible to do this analytically, but it is straightforward to see 
that 

h, J2EI ^ T,Q ^ — (35) 
multiplied by a numerical factor of order unity. In addition we have 

^"^^ ^hli-h ~ -Kp—- ~ "^h^i-h ~ rX. (36) 



d log e d log e 

The leading order term in the denominator of (|33|) is Ji ~ /^"^loge whence 

rl^rl{\oge)-' (37) 



dJEi _ 1^ 
d log e /i Ji 

and our final estimate is 



A^(log6)-^ (38) 



9 loge 

up to a numerical factor of order unity. This dependence on log e will be negligible 
in the regime where e 0. As we discuss in more detail in the next section, this 
provides an important justification for the adiabatic approximation used in this 
paper, which is the context in which we interpret the logarithmic divergences in the 
partition function (||) and the entropy (H). 

We complete this section by commenting on the relationship of this analysis to 
that of 0. We are fortunate in that the terms in H which are dependent on loge 
are very small, leaving in effect 

^ = -7 E {^Eia^i + ^ (1 + e-^'") (39) 

^ E,l,m 

which is to be compared with the expression 

H = -\Y. [(^Eia^Ei + a^EfiEi) E (40) 

^ E,l,m 

which was employed in [0. The only difference between these two expression is the 
e"^'^'' term. As discussed in [0 , the WKB approximation used is only valid for black- 
hole solutions with large radii r/i, for which 5h is very small. The two expressions 
and (|^) are approximately the same in this regime. 



4 Modified Quantum Liouville Equation 



We now turn to the interpretation of the logarithmic divergences in the partition 
function (|^), in the entropy ®, and in the effective Hamiltonian (|38D . As usual, we 
assume that the density matrix for the scalar field is given by the Gibbs formula: 



-PH 



P 



(41) 



where j3 is the inverse of the effective temperature of the EYM black hole given 
in d^). The denominator in (^) ensures probability conservation: Trp = 1 at all 
times. One would expect ( ^11) to be valid (at least approximately), if the black hole 
is (at least approximately) static, as we assume in the adiabatic approximation used 
in this paper. However, we emphasize that we do not have a proof of (^Tj) in this 
particular case. If it is subject to corrections, there may be additional corrections 
to the quantum Liouville equation for p, beyond the one we exhibit shortly. 



If the Gibbs formula (^l]) were exact, we would normally conclude that dp/dt = 0, 
since [p, H] = 0. However, we know from the analysis of Hawking [|T0| that quantum 
effects cause the black hole to decay, changing its effective radius, c.f., (|]): 
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(42) 



The fact that the radius changes as the black hole decays, i.e., changes with time, 
motivates the identification of time t with loge. 



This identification was supported in the two-dimensional string case by explicit 
calculations 0, for example of the black-hole metric @]. In that case, we were able 
to incorporate the effect of back reaction using world-sheet instantons and the Liou- 
ville field, but here we do not have the corresponding formal tools at our disposal. 
Nevertheless, amplitudes and correlation functions in the four-dimensional case ex- 
hibit logarithmic divergences in complete formal analogy with the two-dimensional 
case, supporting the identification of loge with time. 



Since both the Hamiltonian H (|3^ ) and the partition function Z depend on 
loge, this identification implies that, even assuming the Gibbs formula (|41|), there are 
possible sources of time variation in p which are not given by the usual commutator 
term [p, H] : 

dp 
di 



dH dF 



(43) 



where the free energy F is calculated in [|^. Since the normal commutator term 
[p, H] vanishes, the result (|43| ) implies a modification of the quantum Liouville term 
by the addition of a non-commutator term as in (|^), with 



^(1 - 2m',)e-^''^- + —(2-Am'+'^\ e'^^'^ - —rWe~''^ (44) 
1440 ^ e 720 V th ) 2A '^^ ^ ' 

which is our main result. We note that the first term, obtained from the Hamiltonian 
(p!3D, is suppressed relative to the terms obtained from the free energy. Such a 
suppression is necessary for the consistency of our adiabatic approximation. The 
term in (|^) that depends on the energy content of the scalar matter field is the last 
one, which is of 0[/i^], where /x is the mass of the scalar particle. 



The rest of the terms depend on details of the black-hole background, and in 
particular on its hair 0. For extreme macroscopic black holes with rrih » 1 (in 
units of the Planck mass Mp), for instance, the dominant term is the /i^ term. This 
appears not to be the case for non-extreme black holes, where the usual Bekenstein- 
Hawking term seems to be the dominant one, for e ^ 0. However, such terms may 
always be absorbed in a renormalization of the gravitational coupling constants of 
the model 0, as explained at the end of section 2. Therefore, even in this 



case, it is the matter /x^ term in ( PD that determines the order of magnitude of the 
modifications of quantum mechanics p[. 



In general, the quadratic dependence of on the scalar mass fi, divided in order 
of magnitude by just one power of Mp, is the largest that could be expected for 
any such modification of the quantum Liouville equation: a priori, it could have 
been suppressed by one or more additional powers of fi/Mp, or an exponential, or 
even absent all together. We are not in a position to estimate the coefficient of this 
/i^/Mp term. Nor, indeed, can we be sure that such a parametric dependence would 
persist in a complete quantum theory of gravity. However, we would like to remind 
the reader that just such a quadratic dependence also appeared as a possibility in 
a previous string analysis of the quantum Liouville equation. We cannot resist 
pointing out also that such a dependence may not be many orders of magnitude 
from the experimental sensitivity to such a modification of the quantum Liouville 
equation for the K'^ system p, ^ |[. 



Notice also that the main time-dependence in (]4^) comes from the divergences 
occurring when one traces over an infinite number of states in Z. Any explicit time- 



dependence of the Hamiltonian operator H is sub leading at large times (38). This 



is not in contradiction with the adiabatic approximation made above, in which the 
effect of quantum fiuctuations on the Hawking temperature /? was neglected. This 



approximation is valid for macroscopic black holes, for which the WKB approxima- 
tion method we have used to derive expressions for the entropy and the free energy 
is vahd 0. 

It is easy to check that the e-dependence in (^) , combined with the Gibbs formula 
(^TD, reproduces the e-dependence of the entropy found previously in 0. The 
interpretation developed earlier of the e-dependence as a time-dependence implies 
that the entropy of the scalar field increases with time: 

^ = - — (1 - 2m'^)e-^'^ - - — {2-Am'^ + '^ e'^''^ + -r^/x^g-^'^ (45) 
dt 360^ e 180 V th j 12 ^ ' 

Thus the (f) state evolves to become more mixed, as a result of the non-commutator 
term (|4^) in the quantum Liouville equation (0). This can be thought of as being 
due to a change in the entanglement of the external field with the unmeasurable 
modes interior to the black hole. The amount of this entanglement depends on the 
quantum numbers of the EYM background, as can be seen in @. The more hair it 
has, the smaller the e-dependence in and, correspondingly, the slower the rate 
of information loss in (^). 



Although our analysis breaks down in the limit of an extremal black hole, the 
results (?4) and (^) suggest that there may be information loss even in this case. 
This is because it is no longer true in general at the quantum level that the entropy 
is proportional to the area of the horizon (^. There will be information loss (entropy 
increase) even in the absence of any finite-temperature effects, if there is entangle- 
ment with modes beyond the horizon at the quantum level, as we have illustrated. 
This observation is related to the phenomenon of entropy generation during inflation 
in cosmology, which may also be regarded as a non-equilibrium process associated 
with information loss beyond the Hubble horizon. In the context of non-critical 
strings 0, we have discussed in ref. 



ICl how such an information loss can lead to 



a stochastic framework for time evolution, during such non-equilibrium processes. 
The stochasticity of the time evolution in the Liouville string |16], where the time 



variable is identified with a RG scale, can be derived from some specific proper- 
ties of the RG evolution in two-dimensional spaces (world sheets) |T^, |]. One can 
hope that a similar framework may be developed here, identifying a renormalization 
group (UV) scale loge with time. However, we are not yet in a position to prove 
that a similar stochasticity characterizes the RG evolution in the four-dimensional 
case. However, the presence of logarithmic infinities in the (entanglement) entropy 
of black holes does seem to be a generic phenomenon, independent of the dimension- 
ality of space-time, in view of the fact that they are present even in two-dimensionsal 
models 



The above analysis is not complete yet, as a result of the absence of a satisfactory 
treatment of quantum gravity. In the absence of such a treatment, the results of this 



paper can only be regarded as indicative. However, we think that they constitute 
interesting circumstantial evidence in favour of the picture advanced previously ||, 
1^ , namely that microscopic quantum fluctuations in the space-time background may 
induce a loss of quantum coherence in apparently isolated systems. Moreover, the 
magnitude of that we find is consistent with previous string estimates, and may 
not lie many orders of magnitude beyond the reach of particle physics experiments 
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